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ABSTRACT 

A generalized plane-strain elastic -plastic analysis and a creep analysis of a long 
hollow (or solid) cylinder are presented. These analyses use the method of successive 
approximations for solution. The equations are converted into finite difference formu- 
lation and programmed on an IBM 7094 Mod n digital computer. The cylinder is sub- 
jected to thermal gradients and pressure loads. The resulting stresses and strains cal- 
culated are printed as the output. A sample problem is included, along with program 
operation instructions and listings. Examples of components for which the analysis can 
be used are reactor pressure shells, fuel elements, and heat- exchanger tubes with sym- 
metrical temperature gradients and constant pressure loads. 



A DIGITAL COMPUTER PROGRAM FOR DETERMINING THE ELASTIC- 
PLASTIC DEFORMATION AND CREEP STRAINS IN 
CYLINDRICAL RODS, TUBES AND VESSELS 
by Richard L. Puthoff 

SUMMARY 

Many design applications, such as nuclear reactors, often encounter both high heat 
fluxes and internal heat generation. As a result, these structures are subjected to large 
thermal gradients in addition to the pressure loads. These thermal stresses, coupled 
with high temperatures of operation, necessitate that the material operate near its maxi- 
mum load-carrying capacity. Therefore, the materials are often operated into the plas- 
tic range. 

A generalized plane-strain elastic -plastic analysis is presented of a long hollow (or 
solid) cylinder with thermal stresses and constant pressure loads. The analysis uses the 
method of successive approximations for solution and applies a one-step loading path 
called "deformation theory" for the calculations. The equations are converted into finite 
difference formulation and programmed on an IBM 7094 Mod II digital computer. Stresses 
and strains due to the thermal gradients and pressure loads are printed as the output. 
These values represent the elastic -plastic stress state existing on the cylinder when the 
loads are initially applied (time, 0). Further deformation by creep may then be deter- 
mined. 

The generalized plane-strain creep analysis is similar to the elastic -plastic solu- 
tion, except that the equivalent creep strain replaces the equivalent plastic strain and 
plastic flow occurs at all mesh points. Here, however, the loading path is no longer a 
one-step process, and the "incremental theory" is applied to the calculations. These 
equations have also been programmed on an IBM 7094 Mod II digital computer. 

A sample problem of a hollow steel cylinder (appendix A) has been included to demon- 
strate the use of the code. Program operation instructions and listings are included in 
appendixes B and C. 



INTRODUCTION 


In nuclear reactor applications high temperatures of operation, coupled with internal 
heat generation, result in unique design problems. Examples of components with these 
problems are pressure shells, fuel elements, and heat- exchanger tubes. In the case of 
the pressure shell surrounding the core; not only is heat generated within the shell by 
gamma radiation emitted from the core, but the usual radiative and conductive heat 
sources are also present, as in any vessel surrounding a heat source. In the case of the 
fuel elements, high heat generation occurs because of neutron fissioning. The heat is 
then carried away by coolant flow through the core. In both cases large temperature gra- 
dients occur, resulting in high thermal stresses. These thermal stresses, coupled with 
high temperatures of operation, necessitate that the material operate near its maximum 
load-carrying capacity. Therefore, these materials are often operated into the plastic 
range. 

The solution of elastic -plastic flow problems involves the use of equilibrium and com- 
patibility equations combined with nonlinear stress-strain relations instead of the linear 
Hooke's Law. The solution to the resulting equations requires the use of the method of 
successive approximations. This method has been used by Millenson and Manson (ref. 1) 
in connection with the rotating disk and by D, F. Johnson (ref. 2) in analyzing thin-wall 
cylinders and spheres subjected to internal pressure and nuclear radiation heating. 

The method of successive approximations is also applicable to the solution of creep 
problems. This method can use various assumptions with regard to cumulative creep, 
such as the strain- hardening, time -hardening, or life-fraction rules (ref. 3). The time- 
hardening rule was used in this report. 

A digital computer program is presented, for use on an IBM 7094 Mod II computer, 
that calculates stresses and strains in a long hollow (or solid) cylinder undergoing elastic- 
plastic deformation and creep. FORTRAN IV compiler language has been used. All equa- 
tions which were used and converted into finite difference formulation are presented. A 
sample problem and program listings, together with the operating instructions, have been 
included in the appendixes. 


SYMBOLS 


A creep law constant 
a radius, in.; cm 

b radius, in. ; cm 

2 

E modulus of elasticity, psi; N/cm 


2 



AH material constant in creep law 
h radial increment in finite difference, r n - r n _ ^ 

K constant 

m creep law exponent 

N creep law exponent 

n creep law exponent 

2 

P internal pressure, psi; N/cm 

Q/A surface heat flux, Btu/(hr)(ft^); W/cm^ 

R universal gas constant 

r radius, in. ; cm 

T time, hr 

AT time increment, hr 

t temperature, °F 

a coefficient of thermal expansion, in. /in. /°F; cm/cm/K 

Ae incremental strain, in, /in. ; cm/cm 

e rate of strain, in. /in. /hr; cm/cm/hr 

e ec equivalent creep strain, in. /in. ; cm/cm 
Ae ec incremental equivalent creep strain, in. /in. ; cm/cm 
e et equivalent total strain, in. /in. ; cm/cm 
e equivalent plastic strain, in. /in. ; cm/cm 

r 

e r radial strain, in. /in. ; cm/cm 

Ae r incremental radial strain, in. /in. ; cm/cm 

e, longitudinal strain, in. /in. ; cm/cm 
z 

Ae_ incremental longitudinal strain, in. /in. ; cm/cm 
€q tangential strain, in. /in. ; cm/cm 
Ae q incremental tangential strain, in. /in. ; cm/cm 

p Poisson's ratio 

9 

a equivalent stress, psi; N/cm 

V 

n 

cr r radial stress, psi; N/cm 
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o 

(J z longitudinal stress, psi; N/cm 

2 

tangential stress, psi; N/cm 
Subscript: 

n radial mesh points 
Superscripts: 
c creep 

p plastic flow 


ANALYSIS 

The analysis of a hollow (or solid) cylinder (fig. 1) presented in this section consists 
of both an elastic -plastic analysis under the initial loads and a creep analysis. They are 
developed separately and programmed separately. Both analyses use the method of suc- 
cessive approximation for the solution. In the case of elastic-plastic flow, the deforma- 
tion theory is used (ref. 4), whereas in the creep solution the incremental theory is ap- 
plied (ref. 4). 

To develop the analysis the following assumptions were made: 

(1) Axial symmetry exists. 

(2) Steady-state heat flow exists. 


x 3 
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(3) Materials are isotropic. 

(4) The von Mises’ yield criterion is used. * 

(5) Volume change due to plastic flow is zero. 

(6) The cylinder consists of a strain-hardening material. 

(7) Generalized plane strain exists. 

(8) Hie cylinder is sufficiently long that Saint Venant’s principle applies. 


Elastic- Plastic Analysis 

The following is a derivation of the equations for calculating the elastic-plastic 
stresses and strains in a hollow (or solid) cylinder with a symmetrical temperature gra- 
dient and internally or externally pressurized. The equilibrium equation is 



dr r 


The compatibility equation is 



dr r 


and the stress-strain relations are 

€ r = ^ [ a r - + t “ t + 

c 9“^[ tr S- + “z)] +<rt + 

£ z =j[ < ’ z -' 1 <‘ T e + a r>] +at - s ? 

where by assumption (5) the plastic strains are 
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( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 


1 This yield criterion is also referred to as the distortion energy theory. It assumes 
that yielding begins when the distortion energy equals the distortion energy at yield in 
simple tension. 
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( 6 ) 


P . .P + P P - o 
e r + e 6 + e z “ U 


Substituting equations (3) to (5) into the compatibility equation (2) gives 


J i ff 0 ' e " r ‘ ' (<T r + °e> + I* 0 * + >>(-4 - e ?) + + £ 1 


e P _ e P 

1 + a /„ _ * , e r € 0 

= — ( a r - V + 

Er r y r 


( 7 ) 


For the generalized plane-strain problem, since e z is constant, 

d(-/xe z ) 


dr 


= 0 


For the solid cylindrical configuration, 


r-. 


r dr = 0 


and therefore, 


e 


z 




Er dr 


For the case of the hollow cylindrical configuration, 



o r dr = 
z 2 


( 8 ) 


and e is changed accordingly, 
z 

This equation expresses the axial strain as a function of stresses o r and o^, plastic 
strains and e^, temperature t, and material properties E and a. 

The method used for calculating the stress-strain state of the cylinder was the strain- 
strain method because of its strongly convergent characteristic. This method relates the 
equivalent total strain with the equivalent plastic strain. 
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When a material is subjected to the stress state (cr r , <Jq, a z ) the strain components 
can be determined by application of the following basic assumptions: 

(1) The principal plastic shear strains are proportional to the principal shear stresses 



°r “ a 6 



(9) 




= 0 


(3) A universal relation exists between the equivalent stress a and the equivalent 
plastic strain e n where 

ir 


a 


e 





e 


P 




( 10 ) 


(ID 


When the strain-strain method is being used, the elastic and plastic components of strain 
are separated by using the concept of equivalent total strain defined by 


e et 



e z>‘ 


( e fl - 



( 12 ) 


The equivalent plastic strain is related to the equivalent total strain by 


£ et 


+ 2(1 + M ) g e 
3 E 


(13) 


From this relation a uniaxial curve (fig. 2(a)) can be replotted as a curve of e p against 
e et ( fi §- 3(b)). From equations (9) to (13), the individual plastic strains for the strain- 
strain method become 


e 


P = 
r 



(2e r - - € z ) 


(14) 
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(a) Uniaxial stress -strain curve. (b) Strain -strain curve. 

Figure 2. - Equivalent plastic strain and equivalent total strain obtained from uniaxial curve. 


4=jE_(2€ e -€ r -£ z ) (15) 

35 et 

and can be computed from equation (6). These equations express the radial and tan- 
gential plastic strains as a function of the equivalent total strain, the equivalent plastic 
strain, and the total strains. 


Finite Difference Formulation 

In the finite difference method, a number of discrete mesh points are chosen at inter- 
vals along the radius which are not necessarily equal. There are thus n + 1 mesh points, 
the first at the center of the cylinder (or inner radius) and the last at the outer radius. 
Using the method of middle differences gives the typical equations 

da _ °n ~ CT n- 1 
dr r n‘ r n-l 


2 
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Equation (7) can be restated as 


d 

1- M 2 

dr 

E 


(iL + iL. U + (l + n)at + 


(E E, 


Applying equations (16) results in 
2 /.. 2 \ 


i±Jicr r -i±Jicr fl+ 5 e ' 
Er r Er ° r 


’ T K" 


<rp 

e r e 0 


:) 

(17) 


^ 1-lf.. (v Q -On h) - + J£— (a - O 1 ) 

Eh \e E / h r? n r > n- l 

n \ 7 n 

+ (1 + l*)[ (at) n - (at)^] - -t (€ z n - e Zi n .j) 




+ JL Lp . A + £ p )+e P. .Uife -if A L i±Jt (a + „ ) 

h \ ^ n r, n- 1 0 ? n- Xj ^ l 0, n 9 y n- 1 J 2Er ^ n- 1 

n n ^ 

- i-i-H (a„ + A + -L/gP + gP - ~(e? + € § ^ ( 

2Er ■ ? b U“ 1 2 j> 1 b r ? n~ 1 1 ^ 0, n b — ly 


Also applying equations (16) to equation (1) gives 


°r, n ~ °r, n- 1 + a r, n + a r, n- 1 u 0,n ' ^n-l _ 

2r 


r n ' r n-l 


2r 


(19) 


Now, by collecting all the n and n - 1 terms, equations (18) and (19) can be restated as 
follows: 

( 20 ) 
( 21 ) 

where 


C' cr + D* da _ = F' cr 1 + G'a ft „ , + H’ + P’ 
n r, n n 0, n n r, n-1 n 0, n-l n n 


C ct „ - D „ = F ct _ ! + Ga. _ « 
n r, n n 0, n n r, n- 1 0, n- 1 
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C = — + — i— 

h n 2r n 

n n 


C;t - zR - iL _ 1 + M 

n hE n hE n 2E r 
n n n n n n 


D n = 


2r 


n 


F = 
n 


h n 2r 1 
n n- 1 


G n~ 


2r. 


n-1 


D* 

= 1 


1 + /i. 

n 

h E 
n n 

V. 

2Er 
n n 

f’ = 


. * 2 

1 + M 

r n 

i 

1 hE 

1 n n- 

i 2E .r 
-1 n-1 n 

G* = — 

1 

M 2 . 

1 + ii 

n h 

E . 

h E . 

2E .r . 

n n-i 

n n-1 

n-1 n-1 

1 + n 
h 


< at >n-l] 



n 


P’ J.A L+jAjP +1 1 + ^V p +(- 1 - 1 - Mr p , / 1 h\ c P 

11 \ h n 2r n V ’ n \ 2r n h / ^ \\ 2 Vl \/ ^ \ 2 Vl h n) ^ 


y (22) 


Considering the linear nature of equations (20) and (21), it follows that the stresses at any 
mesh point can ultimately be expressed in linear terms of the stresses at any other mesh 
point. Solving for the stresses at the n m mesh point in terms of stresses at the n-1 
mesh point gives 


a r, n *11, n°r, n- 1 + *12, n CT 0, n- 1 + M l, n 


° 9, n *21, n ff r, n- 1 + *22, n CT 0, n- 1 + M 2, n 


(23) 

(24) 


where 
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D' F + D F' 
_ n n n n 

11,11 CD’ + C’D 
xi n n n 


C F’ -C’F 
, _ n n n n 

21,11 CD’ + C’D 
n n n n 


1,n c „ D ; +c ; D n 


D’ G 

+ 

D G’ 

n n 


n n 

c„d’ 

+ 

C’D 

n n 


n n 

C G’ 


C’ G 

n n 


n n 

CD’ 

+ 

C’ D 

n n 


n n 

c n (ir • 

n\ n 

^ P n> 

C D’ 
n n 

+ 

C ; D n 


In matrix notation 

a - La < + M„ (26) 

n n n- 1 n v ' 

where a n , L n , <? n _j, and M n are the indicated matrices. Equation (26) represents a 
linear recurrence relation between the stresses at the n mesh point and the n - 1 mesh 
point by which they can be linearly related to the stresses at the first mesh point. 

ff n - Vl + B n < 27 > 

where A n and B n are yet unknown. Substituting equation (27) into equation (26) gives 

(A - A 1 L n )a 1 = L B , - B n + M n (28) 

n n- 1 n' 1 n n- 1 n n v ' 

Now a^ depends on the boundary conditions and therefore is completely arbitrary; how- 
ever, equation (28) must be true for all values of a^. Therefore, both sides of the equa- 
tion must vanish identically, as the A’s and B’s are independent of the boundary conditions 
and are functions only of the geometry and material properties. Therefore, 

A = A .L (29) 

n n- 1 n 


B — LB ■< + M 
n n n- 1 n 


Now, for the second mesh point, equation (27) gives 


°2 = A 2 ct 1 + B 2 
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and equation (26) gives 


a 2 = L 2 a l + M 2 (32) 

Therefore, Ag = Lg, Bg = Mg, and by the recurrence relations (eqs. (29) and (30)) all 
other A n and B r can be computed. For the solid cylinder, the following boundary con- 
ditions were used: 


ct = -P at r = a 
r, a 


(Ti = o a 1 = <j « at r = 0 
1 9, 1 r, 1 


(33) 


Substituting equations (33) into equation (27) yields 


°9, n “ A 0,n °9, 1 + B 0,n 


CT — A (7^ h + B 

xyn r ,n 0,1 r, n 


Therefore, A r j ^ = 1 and B r ^ = B 0 = 0. 


-B + a 
- r,a r,a 

a M 7 

A r,a 


(34) 


The tangential stress at the inner mesh point is a function only of the coefficients a , 
B„ _ and the radial stress at r = a. After A„ and B are derived from equations (29) 
and (30) and Oq ^ from equation (34), the stresses at every mesh point can be computed 
from equation (27). 


Creep Analysis 

The creep analysis is treated in a manner essentially similar to that of the elastic- 
plastic flow. However, since the loading path is now no longer a one-step process but 
rather an incremental time- loading process, the equations derived in the previous section 
must be altered. The following, therefore, is the analysis of a hollow (or solid) cylinder 
with a symmetrical temperature, internally or externally pressurized and undergoing a 
creep over a predetermined period of time. 
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Modification of elastic-plastic equations . - The equilibrium and compatibility equa- 
tions (eqs. (1) and (2)) remain the same, but the stress-strain relations now reflect the 
incremental creep strain and the total creep strain. These relations are 


1 

E 


E' 


+ at + + Ae£ 

r r 

(35) 

+ at + € ^ + A 

(36) 

c e A c A c 
€ r ~ e G “ A€ r “ Ae 0 

(37) 


where again the incompressibility assumption applies to both the incremental and total 
creep strains. The compatibility equation in terms of stresses becomes 


d 

[\ 2 
1 “ M 

dr 

. ^ 


Hf^ 9 '(f + ^ +(1+A) “ t+ K' 


c c A c A c\ c a e 
e r " e e ” Ae r ' Ae e) + € e + A€ e 


c . c c . c 
e„ + Ae_ - Sq - A €q 


= l±ii (CT - On) +-£ l 


Er 


(38) 


The equivalent stress is calculated by using equation (10). The strain increments due to 
creep are then calculated by 


(39) 

(40) 

Ac ^ = -Ae£ - Ae c $ (41) 

where Ae is obtained from the power law relation between the equivalent stress, 
equivalent strain, and total time T. This will be discussed further in the next section. 

In the finite difference formulation, the middle difference equation (eq. (16)) is again 
applied to equation (38). The geometric relations are then derived. These equations re- 
main identical to equations (22) with the exception of which becomes 


4£ r - ~ < 2a r - ”6 ‘ 

A£ 9 “ ^ < 2<r e - ’ °z> 
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Stress-strain - time relations . - In creep the strain path becomes time dependent. 
The relation between the strain and time can be expressed by three basic laws (ref. 5), 
the time-hardening rule, the strain-hardening rule, or the life fraction rule. 

There are many stress-strain - time relations in existence (ref. 5); however, the 
ones used are often dependent upon the data available. A typical example is the constant- 
temperature stress-strain - time relation. This is one in which the logarithm of the 
linear creep rate against the logarithm of the stress is linear. 

log e = log A + N log a 


Ae ec - Act N A T 

where the values of N and A are based upon data available. Since this relation holds 
only for a constant temperature, similar relations would be necessary for a case where 
a wide range of temperatures are encountered. Also this relation is valid only in the 
secondary creep regime. 

Other creep functions are a power-law relation of equivalent stress, equivalent 
strain, and the total time T during which the stress persists 

and a creep function in which the creep curves at a given stress but different temperatures 
can be brought into coincidence 

<ec = ^ 

where k' = Ke -n This relation was first introduced by Dorn (ref. 6). 


or 


and 
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METHOD OF CALCULATION 


The programming of the elastic-plastic and creep analyses was performed separately. 
The elastic -plastic analysis is used first and represents the cylinder at a time equal to 
zero. Since the temperature increase is not uniform, the volumetric increase is also not 
uniform, which results in thermal stresses. Superimposed on these thermal stresses 
are stresses as a result of pressure loads. The stress state resulting from this initial 
loading and occuring at a time equal to zero then becomes the input to the creep analysis 
for calculating the stress state after a period of time. 

The following is a brief summary of the program technique of calculating both elastic- 
plastic deformation and creep. Program listings are given in appendix C. 


Summary of Method of Calculation for Elastic- Plastic Deformation 

The procedure for computation is as follows: 

(1) Set total plastic strains equal to zero, as starting condition. 

(2) Or set total plastic strains equal to plastic strains of previous iteration. 

(3) Solve for stresses and total strains from the sets of equilibrium equations, com- 
patibility equations, and stress-strain relations, as in any elastic problem. A first ap- 
proximation is thus obtained for the stresses and total strains. This includes the geo- 
metric calculations of equation (22), the matrices L n and M n of equation (25) and A n 
and B n of equations (29) and (30), and finally the stresses using equation (27) with the 
aid of equation (34). 

(4) Calculate the equivalent total strain from equation (12). 

(5) Calculate the equivalent plastic strain from equation (13). 

(6) If 6p is positive, yielding has occurred. Compute plastic strains. 

(7) Compute plastic strain components from equations (14) and (15). Return to step 

( 2 ). 

(8) Repeat process until convergence is achieved. 

(Convergence is required of the total strains e r and eg at each mesh point.) A sample 
problem utilizing this method of calculating elastic -plastic deformation is included in ap- 
pendix A. 


Summary of Method of Calculation for Creep Deformation 


The final calculations of the elastic-plastic solution (which represent the stress- 
strain state at time zero) become the input to the creep calculations. The method then 
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proceeds as follows: 

(1) Set the initial total strains and of the creep analysis equal to the total strains 
and P^ of the elastic-plastic calculations. 

(2) Choose a time increment. (The size of this interval depends on the particular 
problem. ) 

(3) Solve for the stresses and total strains from the sets of equilibrium equations, 
compatibility equations, and stress-strain relations, as in the elastic -plastic calculations. 

(4) Calculate the equivalent stress using equation (10). 

(5) Calculate the equivalent creep strain from equation (43). 

(6) Calculate creep strain increments using equations (39) to (41). 

(7) Repeat steps 3 to 6 until two successive calculations of the plastic strain incre- 
ments Ae£, Ae g at each mesh point remain within the convergence criteria. 

(8) Add incremental strains (eqs. (39) to (41)) to total strain and set incremental 
strains equal to zero. 

(9) Return to step 3 for next time increment. 

(10) The problem is completed when total time has been reached. 

A sample problem utilizing this method of calculating creep deformation is included 
in appendix A. 


CONCLUDING REMARKS 

A computer program for the stress analysis of cylindrical rods, tubes, and vessels 
is presented. The primary advantage of this analysis is that it allows the material used 
in a design application to operate in the plastic range. Other advantages are as follows: 

(1) It takes into account the material property changes across the cylinder wall. 

(2) It calculates stresses and strains in cylinders of composite material. 

(3) Both hollow and solid cylindrical geometries can be analyzed. 

(4) Computer running times are short (less than 1 minute for the elastic-plastic anal 
ysis and 2 minutes for the creep analysis of the example problem). 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 1, 1968, 

126-15-01-04-22. 


16 



APPENDIX A 


EXAMPLE SOLUTION 

The operating instructions for using the elastic-plastic program and the creep pro- 
gram are presented in appendix B. To further clarify the method of using these programs, 
an example calculation is presented. The problem involves a typical high-alloy steel tub- 
ing 2. 50 inches in outside diameter and 2. 25 inches in inside diameter with a temperature 
gradient across the wall. The heat flux Q/A is from the inside wall outward. No pres- 
sure loads exist, leaving only thermal stresses. The following steps were taken in the 
solution of this example problem in the elastic-plastic program: 

(1) Obtain uniaxial stress-strain curves for the material at different temperatures 
(fig. 3). 

(2) Convert each stress- strain curve to an equivalent plastic strain - equivalent total 
strain curve by using equation (13). (See fig. 4 for a typical strain-strain curve. ) The 
stress-strain curve has been approximated by two straight lines. The equations of these 
straight lines in the plastic region form a part of the subroutine PLSTR. 

(3) Divide the 0.125-inch wall into any desired number of mesh points, recording their 
radial location. 



Figure 3. - Stress as function of strain for high- 
alloy steel tubing. 



Figure 4. - Equivalent total strain as function 
of equivalent plastic strain at 1000° F for 
high -alloy steel tubing. 
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Figure 5. - Temperature profile through wall of tubing. 



(4) Using the thermal temperature gradient of figure 5, record the respective tem- 
peratures of each mesh point. 

(5) Complete the program input by recording the values of modulus of elasticity and 
coefficient of thermal expansion at each mesh point for that respective temperature. 

The input is complete, and the program output will be as outlined in appendix B. The 
input- output of the sample problem is as follows: 


NO Of ITERATIONS = 

1 

NO OF STATIONS = 

14 


STATION NO 

CYL RADIUS 

TEMPERATURE 

ALPHA 

6 

1 

1.125000000 

1200.000000 

0.10200DOQOE-04 

22000000.00 

2 

1.135000005 

1100.000000 

0. 10 .100300 16-04 

22500000.00 

3 

1.144999996 

1000.000000 

0.1000000006-04 

2 30000 00. 00 

A 

1.155000001 

910.0000000 

0 • 99000000 IE— 05 

23500000.00 

5 

1.165000007 

812.0000000 

0.98 000000 4E-0 5 

24000000.00 

6 

1.174999997 

720.0000000 

0 • 96 9999 99 6E— 05 

24500000.00 

7 

1.185000002 

625.0000000 

0. 960000000E— 05 

25000000.00 

8 

1.194999993 

530.0000000 

0. 950000 003E-0 5 

25500000.00 

9 

1.204999998 

432.0000000 

0 • 9399 99 99 5E— 0 5 

260000on.no 

10 

1.215000004 

335.0000000 

0 .929999 999E— 0 5 

26500000.00 

11 

1.224999994 

240.0000000 

0 .9299999996—05 

27000000.00 

12 

1.234999999 

145.0000000 

0.9200000026-05 

27500000.00 

13 

1.245000005 

50.00000000 

0.920000 00 2E-05 

28000000.00 

14 

1.250000000 

0 

0 • 920000 00 2E-05 

28500000.00 

POISSONS RATIO = 

0.30000 RAD R 

STRESS = 0 

RAD I STRESS 

= 0 
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STATION 

NO 

TOTAL STRAIN R 

TOTAL STRAIN T 

TOTAL STRAIN Z 


1 


0. 1 802 6044 8E— 01 

0. 5477302186-02 

0.5501926186-02 


2 


0.1582725956-01 

0.5578224436-02 

0. 5501925956-02 


3 


0. 136874190E-01 

0.56 58435646—02 

0. 5501925836-02 


4 


0. 1 1 7882340E-01 

0. 5719 7 6 74 2E -0 2 

0. 550192 58 3E-02 


5 


0.97951 53 056-02 

0. 5763341966—02 

0.5501925716-02 


6 


0. 796385459E— 02 

0.57898967 IE-02 

0. 550192 577E-02 


7 


0.612923561E-02 

0. 5800534516-02 

0.5501925836-02 


8 


0. 43451 1906E-02 

0.5795 85 1476-02 

0.5501925656-02 


9 


0* 255839454E-02 

0.5776428446-02 

0.5501925716-02 


10 


0. 839047774E— 03 

0. 574289 513 E-02 

0.5501925656-02 


11 


-0. 75308233 5E— 03 

0. 5696390176—02 

0.5501925656-02 


12 


-0. 236056492E-02 

0. 563768391E— 02 

0.5501925606-02 


13 


-0.3 91 24448 7E— 02 

0. 556723 110E-02 

0.5501925486-02 


14 


-0.4726 9765 IE— 02 

0.5527686276-02 

0. 5501925656-02 


STATION 

NO 

PLASTIC STRAIN R 

PLASTIC STRAIN T 

CYL RADIUS 


1 


0. 5490465 036-02 

-0.2753320766—02 

1.131161943 


2 


0.4021009146-02 

-0.198813717E-02 

1.141331270 


3 


0.259709 92 9E-02 

-0.1 260947006-02 

1.151478887 


4 


0. 1 72799350E-02 

-0. 8 18294313E-03 

1.161606327 


5 


0. 491611 441E -03 

-0.222649 9336-03 

1. 171714291 


6 


0 

-0 

1.181803107 


7 


0 

-0 

1.191873625 


8 


-0 

0 

1. 201926023 


9 


-0 

0 

1.211960584 


10 


-0. 9507882426-03 

0.5113172166-03 

1.221977606 


11 


- 0. 1 923 32 81 3E— 02 

Q. 1005 823976-02 

1.231978059 


12 


-0. 2902332 876-02 

0. 1488429766-02 

1.241962522 


13 


—0. 384348340E-02 

0. 1941668686-02 

1.251931190 


14 


-0. 433673261E-02 

0.21 76547286—02 

1 .256909 594 


STAT ION 

NO 

RADIAL STRESS 

TANGENT STRESS 

LONG STRESS 

EQUIV STRESS 

1 


-0 

-212363.3438 

-211946.6289 

212155.8242 

2 


-1725.279205 

-179112.4277 

-180432.9785 

178051.5410 

3 


-3127.030182 

-145178.2773 

-147947.2969 

143456.1582 

4 


-4222.230225 

-113921.4355 

-117859.3408 

111720.5000 

5 


-5016.575806 

-79450.01172 

-84276.15820 

76960.27734 

6 


-5510.375732 

-46481.12109 

-51908. 26563 

43936.54150 

7 


-5711.390991 

-12032.56641 

-17775.04102 

10451.45972 

8 


-5618.997131 

22837.67773 

17072.20996 

26056.86670 

9 


-5232.874634 

59127.80469 

53637.74902 

61798.97363 

10 


-4554.446777 

95408.59766 

90496.52734 

97600.00293 

11 


-3596.340088 

130354.2441 

126315.3652 

131977.8320 

12 


-2364.607666 

166829.1211 

163957.3105 

167776.6777 

13 


-859. 6215820 

203318.0195 

20 1911.4375 

203478. 5C39 

14 


-0 

224813.7637 

224249.0117 

224532.4824 

1.000000000 

-1.000000000 

-1.000000000 

-0.5501925776-02 


STATION 

NO 

ELASTIC STRAIN 

PLASTIC STRAIN 

P PRIME 


1 


0. 835761 137E-02 

0.54 90 45 90 9E —0 2 

500000000000 


2 


0. 685824704E— 02 

0*4021082026 —02 

-0.91 3353693E-01 


3 


0. 5405567126-02 

0.2 597455696—02 

-0.8928876006-01 


4 


0.41201678 7E-02 

0.1728794776-02 

-0.5427 180426-01 


5 


0.277910713E-02 

0.4923367816-03 

-0.77 37768636—0 1 


6 


0. 1 55420994E-02 

0 

-0.3002731716-01 


7 


0.36231 543 7E-03 

0 

0 


8 


0.8855883936-03 

0 

-0 


9 


0.2 0599553 8E-02 

0 

-0 


10 


0.3191 933716-02 

0.951690250E-03 

-0.6491750576-01 


11 


0. 4236304146—02 

0.1923999116-02 

-0.6558899126-01 


12 


0. 52 87480726-02 

0.2902644516-02 

-0.6612571986-01 


13 


0. 6298112396— 02 

0* 3843542606—02 

—0.6406 21968E— 0 1 


14 


0. 6 82 7853916— 0 2 

0.4336731976-02 

-0 .67406682 3E-01 
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STATION NO Ml M2 B1 B2 


1 

500000000000 

500000000000 

-0 

0 

2 

158.6089420 

36162.81934 

158.6089420 

3616 2.81934 

3 

157.8848305 

36313.54590 

633.04 19922 

72956.88867 

4 

142. 7867050 

33126.49805 

1406.085037 

107023.4395 

5 

153.3958225 

35894.60303 

2471.819702 

144256.0762 

6 

143.7804546 

33932.21973 

3829.736908 

179939.3086 

7 

147.0413113 

34995.81299 

5471.903687 

217056.8770 

8 

145. 8653069 

35007.70654 

7398.780457 

254552.2969 

9 

148.9077873 

36035.66455 

9610.550659 

293425.1 719 

10 

146. 0668392 

35640.28906 

12105.66260 

332247.6367 

11 

137. 9669018 

33939.89014 

14871.3769 5 

369695.1914 

12 

141.6810799 

35136.88916 

1790 1.54907 

408633.4414 

13 

139.2829418 

34820.71631 

21195. 72046 

447548.3750 

14 

37.30787516 

18691.26318 

22950.60693 

472329.9336 

NO OF ITERATIONS = 

500 

NO OF STATIONS ■* 

14 


STATION NO 

CYL RADIUS 

TEMPERATURE 

ALPHA 

E 

1 

1.125000000 

1200.000000 

0 • 10 2000 00 OE- 04 

2 2000000,00 

2 

1.135000005 

1100.000000 

0.10100000 IE- 04 

2 2500000. On 

3 

1.144999996 

1000.000000 

0.1300000006-04 

23000000.00 

4 

1.155009001 

910.0000000 

0 .990000 00 IE-05 

23500OO0.00 

5 

1.165000007 

812.0000000 

0 • 98 0000 00 4E- 05 

2400000O.00 

6 

1.174999997 

720.0000000 

0. 969999 99 6E— 05 

245000 00, 00 

7 

1.185000002 

625.0000000 

0. 960000 000 E-0 5 

25GCOnoO.OO 

8 

1.194999993 

530.0000000 

0. 950000003E— 0 5 

25500000.00 

9 

1.204999998 

432.0000000 

0. 93999999 5E-0 5 

26000000.00 

10 

1.215000004 

335.0000000 

G. 9299999996-05 

26500000.00 

11 

1.224999994 

240.0000000 

0 • 92 9 9 99 9 9 9E— 0 5 

27000000.00 

12 

1.234999999 

145.0000000 

0.92000000 2E-05 

27500000.00 

13 

1.245000005 

50.00000000 

0. 92000000 2E-0 5 

28000000.00 

14 

1.250000000 

0 

0.92000000 2E-0 5 

28500000.00 

POISSONS RATIO = 

0.30000 RAD R 

STRESS = 0 

RAD I STRESS 

= 0 


STATION NO 

TOTAL STRAIN R 

TOTAL STRAIN T 

TOTAL STRAIN l 

1 

0.2 32042 02 IE— 01 

0.53030959 IE -02 

0. 550127635E-02 

2 

Q. 1 96806409E— 01 

0* 5445 36702E-02 

0.550127 64 6 E-0 2 

3 

0. 1 62 504 72 5E— 01 

0.555478351E-02 

0.550127 629E-02 

4 

0* 1 34 996 61 9E— 01 

0. 56 3553 3 60E— 02 

0. 550127 635E-02 

5 

0.1 041 54787E— 01 

0.568985 9226-02 

0. 550 127 62 9E- 02 

6 

0. 8071 362 04E -02 

0. 5720 146 06E -02 

0.550127 62 9E-02 

7 

0* 6233 72 1 80E -02 

0.573226705E-02 

0.550127 629E-02 

8 

0. 444671 01 7E -02 

0.572901 760E-02 

0.550 12762 3E- 02 

9 

0. 26572 096 IE -^02 

0. 571098 103E-02 

0. 55012761 7E-02 

10 

0.163131 190E— 03 

0. 567562436E— 02 

0.55012762 3E-02 

11 

-0. 23000548 OE -02 

0. 562060997E-0 2 

0.550127 623E-02 

12 

-0* 4787071 5 5E -02 

0. 55464445 6E -02 

0 .5501 27 629E-0 2 

13 

-0, 71 8751492E— 02 

0.5453 83984E-02 

0 .550 12762 3E-0 2 

14 

-0. 84474001 8E-02 

0. 5400 75 92 8 E-0 2 

0.55012761 7E— 02 

STATION NO 

PLASTIC STRAIN R 

PLASTIC STRAIN T 

CYL RADIUS 

1 

0. 893076 870E -02 

-0* 4539 95 05 2E -02 

1.130965 963 

2 

0. 65 82 1238 7E -02 

-0. 331048 87 LE— 02 

1.141180485 

3 

0. 430534745E— 02 

—0. 2 1 365 602 7E— 02 

1.151360214 

4 

0. 280529633E-02 

-0. 1 36703295E— 02 

1.161509037 

5 

0. 893820659E— 03 

—0. 420681 74 3E— 03 

1.171628669 

6 

0 

—0 

1.181721151 

7 

0 

-0 

1.191792712 

8 

-0 

0 

1.20 1846153 

9 

-0 

0 

1.211881712 

10 

-0. 13 48091 4 IE -02 

0.7065374086-03 

1.221895874 

11 

-0. 285922 867E-02 

0.1462 16770 E-0 2 

1.231885225 

12 

— 0.43 796611 IE— 02 

0* 2204219 75E -02 

1.241849840 

13 

-0. 5 840793 95 E— 02 

0* 2903 9 89 72 E— 02 

1.251790017 

14 

— 0. 660631 910E— 02 

0. 326732 561 E— 02 

1.256750926 


20 



STATION NO 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

24*00000000 

STATION NO 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

STATION NO 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


RADIAL STRESS 

TANGENT STRESS 

LONG STRESS 

EQUIV STRESS 

-0 

-74979.28516 

-74149.11328 

74567.85156 

-664. 364830 Q 

-75829.33203 

-75534.03809 

75017.94336 

-1324.145630 

-76586.88086 

-76963.32520 

75451.84961 

-1941.719864 

-68681.23438 

-69820.54883 

67316.57031 

-2500.925873 

-65478.51660 

-67991.67969 

64271.19482 

-2959.774567 

-47271.15381 

-51396.00781 

46 511.30518 

-3189.023010 

-12832.38379 

-17274.51367 

12472.56042 

—3124.469849 

2 2 02 8.4 8438 

17561.25000 

23243.61938 

-2765.802094 

58309.62695 

54115.53027 

59090.26709 

-2222.319672 

68261.56055 

66032. 24902 

69396.25781 

-1627.158264 

73123.53223 

71997.25781 

74194.1503 9 

-1000.936218 

78314.92383 

77968.15430 

79143.24219 

-342.7797241 

83580.01660 

83894.84766 

84080.86426 

-0 

87130.21387 

87762.55664 

8744 8. 31836 

0. 14249762 0E-04 

0. 8 26648 12 5E -05 

-0.550 127 62 9E- 02 


ELASTIC STRAIN 

PLASTIC STRAIN 

P PRIME 


0.11 8685 322 E— 01 

0. 8931 16 137E-02 

500000000000 


0. 947157713E— 02 

0.6582145526-0 2 

-0.146176621 


0.71 48343 78E-02 

0. 430537 6 85 E -02 

-0.143307388 


0. 528805936E— 02 

0. 2805 590 69E— 02 

-0 • 942491 25 6E— 01 


0.32151 0996E— 02 

0. 8943 31373E-03 

-0.121218413 


0.1 645289558-02 

0 

-0 .55 69 8230 9E- 01 


0. 43237991 5E-03 

0 

0 


0. 789975224E-03 

0 

-0 


0. 196966564E— 02 

0 

-0 


0. 361 82 706 8E— 02 

0. 134860998E— 02 

-0.907458374E-01 


0. 52 41 1 04 8 7E -02 

0. 285946857E-02 

-0.100837685 


0.687398674E-02 

0.4379681606-02 

-0.101985931 


0. 844340504E— 02 

0.58 408 1002E— 02 

—0.9899 53 279E— 0 1 


0. 9265 7702 OE— 02 

0. 6606 43 203 E -02 

-0.104263035 


Ml 

M2 

B1 

B2 

500000000000 

500000000000 

-0 

0 

0.780969054 

178.0608482 

0.780969054 

178.0608482 

1.08 6252 972 

249.8384209 

3.425862819 

430.2653008 

38.32816172 

8892.128784 

45.47364235 

9327.938232 

17.3 5 4 5 83 02 

4060.970306 

143.0075321 

13505.55957 

80.51584 721 

19001.75781 

337.9459915 

32671.27661 

147.0413113 

34995.81299 

759.4712524 

68052.39746 

145. 8653069 

35007.70654 

1471.726669 

103840.1387 

148. 9077873 

36035.66455 

2474.977203 

141033.1738 

38.19695711 

9320.052 490 

3659.877716 

151882.4980 

16.83815813 

4142.190857 

4893.252502 

157627.8203 

17. 90616775 

4440.727112 

6154.450928 

163688.9570 

17.92784357 

4481.958435 

7444.316467 

169810.6113 

4. 785164595 

2397.369751 

8103.187988 

174520.9336 


Both the elastic solution (first iteration) and the plastic solution (convergence) after 
24 iterations are shown. The plastic strains indicated in the elastic solution must be 
ignored, as they represent only the initial plastic calculations and the elastic solution 
was not yet adjusted. Although the sample problem presented is that of a tube containing 
a single material, the application can be extended to composite materials (ref. 7). It 
may be desirable to subscript the Poisson’s ratio variable if large differences exist be- 
tween the composite materials. Also where two dissimilar materials have large differ- 
ences in properties, it may be necessary to select finer mesh points at these inner faces 
to improve convergence of the solution. 

The following steps were taken in the solution of the example problem in the creep 
program: 
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(1) Obtain creep equations similar to equation (43). These equations form a part of 
the subroutine PLSTR. 

(2) Complete the input. The creep input is similar to that of the elastic-plastic pro- 
gram except for the plastic strains and the variable PV. These values are obtained di- 
rectly from the output of the elastic-plastic solution. 

With the input complete the output will be as outlined in appendix B. The input-output 
of the creep portion of the sample problem at the end of 692. 5 hours is as follows: 


NO 0 6 l T6K AT IONS “ 

10 

NO OF STATIONS = 

14 



STATION Nil 

CYL RADIUS 

TEMPERATURE 


ALPHA 

E 

1 

1. 125000000 

1200.000000 


3 .1320030 00 E-04 

22000000.00 

? 

i. 135000005 

1103.003000 


3.101003001E-04 

22500000.00 

3 

1.144999996 

L003. 000030 


0 .100003000E-04 

23000000.00 

4 

i. 155000001 

910.0003000 


0.990003001 E- 05 

23500000.00 

5 

1.165000007 

812.0003000 


0.9800030046-05 

24000000.00 

6 

1. 174999997 

720.0000030 


3 .9699999966-05 

24500000.00 

7 

1.185000002 

625.0000000 


0.9630 3 3300E-05 

25000000 .00 

8 

1.194999993 

530.0000000 


0.950 00 0003E-05 

25500000.00 

9 

1. 204999998 

432.0003000 


3 .9399999956-05 

26000000.00 

10 

1.215000004 

335.3300003 


0 .9299999 99 E-05 

26500000.00 

l 1 

1.224999994 

240.0033030 


0.929999999E-05 

27000000.00 

12 

1.234999999 

145.0000000 


0.9230030026-05 

27500000 .00 

I 3 

1. 245000005 

50.00000000 


0.92 3 33)30 2 E-05 

28000000.00 

14 

1.25 000 00C0 

0 


0.923003302 E-05 

28500000.00 

STATION NO 

PLASTIC STRAIN R 

PLASTIC STRAIN T 

P 

PRIME 


l 

0 • 893076998E -02 

-0 • 45 3 989999 E -0 2 


3 


? 

0.65 62 120 03 E-02 

-0. 331353000E-02 


-3.146176599 


3 

0.430535001E-02 

-0. 2136500006-02 


-3.143307395 


4 

0.2805289995-02 

-0. 1367030 03 E-02 


-0.9424910 04E— 0 1 


5 

0. 8938199975-03 

-0.420699998E-03 


-0.121218398 


6 

0 

0 


-3.556982001E-01 


7 

0 

0 


0 


8 

0 

0 


0 


9 

0 

0 


0 


10 

-0. 1348100015-02 

0. 706539999E— 03 


-0.937458002E-01 


l 1 

-0.2859200005-02 

0. 146 2 1 70 00 E— 0 2 


-0.133837700 


12 

-0. 4379700005-02 

0.2204220016-02 


—0 .101985 894 


13 

-0.5840800005-02 

0. 290398000E— 0 2 


-0.989953000E-01 


14 

- 0. 66063 00 Cl 5 — 02 

0. 326731999E— 02 


-0.134259998 


POISSONS RAT I rj = 

0.30000 KAO 

K STRESS = 0 


RAD I STRESS 

= 0 


STATION NO TOTAL STRAIN R TOTAL STRAIN T TOFAL STRAIN l 


1 

0.2 69003077E-U1 

2 

0. 2333164 96E-01 

3 

0. 198579582E-0L 

4 

0. 15 1487 52 4E -01 

5 

0. 1 19721447E-0L 

6 

0. 8t>42 42 3 78E - 02 

7 

0.b785782b9E-02 

8 

0. 49 875 56 02 E— 02 

9 

0. 3L8581 9C8E-02 

10 

0.679689656E-03 

I 1 

-0. 1794562796-02 

12 

- U. 42 92 900 97 E-02 

13 

-0.67 04851 52 E -02 

14 

-0. 79 7102 71 7E — 02 


0. 45 143305 5 E-02 3 . 49 25 83 3 14E-02 
0# 46993 352 0 E-02 3. 49239155 IE-02 
0.4849444 3 3E -02 3 .49 2599356E-02 
0. 496? 73 9 7 IE-02 3 .49 18 52 48 2JE-0? 
0.5337089466-02 3 . 49 1 899 9 8 5 E-02 
0.53 78099496-02 0 . 49 3273 5646-02 
0.51 00 3 7 085 E-02 3 . 493303 403E-02 
0. 5 L 0698 03 9 E-02 3 . 493 302 3666-02 
0.5 098 5 4 54 3 E-02 3 . 49 3305 258E-02 
0.507 250 120 E-02 3 .493303 896E-02 
0 . 50 2 692 76 3 E-0 2 3 . 49 33 11 45 IE- 02 
0. 4962 1234 IE -02 3 . 49 33 14693 E-02 
0. 4878982 236-02 0 .4933186228-02 
0.483069 52 8 €-02 0 .493321 562E-02 
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ML RAO IUS 

1*130083382 
1 . 140338764 
1.150557488 
1.160736635 
1.170872688 
1.130971146 
1.191048294 
1. 20 1107159 
L.21 1 148039 
1.221167475 
1.231162205 
1.241132423 
1.251078501 
1.256042525 


ST AT I i JIM NO 

PLASTIC STRAIN R 

PLASTIC STRAIN T 

l 

0. 1435952665-01 

-0 • 73 65 5 5 38 7 E— 02 

? 

0. 1193191195-01 

-0.606 22 2 82 3 E -02 

3 

0.95 78667 705-02 

-0.481 3548 76 E-02 

4 

0. 4639 92 445 E -02 

-0.22862 9922 £-02 

5 

0.2 582 0882 3 E-02 

-0. 124478193E-02 

6 

0. 1 463 608 Cl E -04 

-0. 6 80 6 2 543 6 E -05 

7 

0. Ill 59 805 8E -05 

-0 • 464 183906E-06 

fi 

- 0. 2146322305-07 

0 . 1 5 815 8 00 5 £-07 

9 

-0.3 58 742 3526 -05 

0. 201681 35 4 E -05 

10 

-0. 1355526795-02 

0. 710432992E-03 

1 1 

-0.2 86912 062 E — 02 

0. 146722507E-02 

12 

- 0. 439? 75 993 E -02 

0.2210 7645 9 E-02 

13 

-0. 585764 14 OE -02 

0. 2912 32290 E-02 

14 

-0.6625996265-02 

0. 32 7 70 804 2E-02 


J ION MCI 

I NCR PLAS SIR R 

INCR I 

1 

0. 1588750825-04 

-0, 

2 

0. 1 549 9 7 94 7E -04 

-0, 

3 

0. 1512 969985-0^ 

-0, 

4 

0.2932070545-04 

-0, 

5 

0.2771683 67c -04 

-0. 

6 

0.4587371485-06 

-3 

7 

0. 4343 72107c -07 

-0 

8 

-0. 2421 6 54 70c - L 1 

0, 

9 

-0. 5526859735-07 

0, 

10 

-0. 12681 06C6E-06 

0 

1 l 

-0. 1 74871 8665-06 

0 

12 

-0. 2364972345-05 

0 

1 3 

-0.3117318355-06 

0 

14 

- 0. 3681 06 01 95 -06 

0 


STATION NO 

RAOUL STRESS 

TANGENT STRESS 

LUNG STRESS 


1 

-0 

-10774.02014 

-10 1 18.16797 


2 

- 96. 69426823 

-10728.58154 

-10 16 1.55469 


3 

- 190. 9273396 

-1 0683.74597 

-10205.00293 


4 

-492.2262650 

-58021.31055 

-58209. 18945 


5 

-991.0944977 

-57608.57422 

-580 74.25879 


6 

-1520.668396 

-68277.66211 

-70996.45996 


7 

-1942.465012 

-34322.64941 

-37537.51563 


8 

-2069.390594 

272.4332275 

-3 139.496094 


9 

-1901. 501617 

36239.20801 

32938.42676 


10 

-1549. 161514 

45880.91797 

44369.05469 


11 

-1145.579605 

50463.48877 

49871.96875 


12 

-711.3513031 

55370.56738 

55366.18457 


l 3 

-245. 7683792 

60350. 62012 

60807.34375 


14 

-0 

63536.29199 

64206.68066 


4.000 000000 

0. 2412 83 3 l 7E -05 

0 . 3647 72 49 8 E -05 

-0 .493302877E-02 

692.4847336 


STATION NO 

EOUiV. STRESS 

PLASTIC STRAIN 

P PRIME 


1 

10461.55 042 

0. 159 1 1 01 6 I E -04 

0 


2 

10360. 04419 

0. 155172 7 47 E -0 4 

-0.142484576 


3 

10261. 85168 

0. 15142 1016E -04 

-0.140488945 


4 

5 7924.26 807 

0 • 293 2 1 3 3 5 9 E -0 4 

-0.309201114 


5 

56851.89453 

0. 277176 03 9 E -04 

-0 .1283050 25 


6 

68157.24414 

0 . 459 3 1 2 27 OE-06 

-0.162168406 


7 

34101. 54639 

0. 435828 163E-37 

-0 *8 59 0 38 708E-03 


8 

3022.467316 

0. U5109960E-10 

-0 .69 2464 42 8E— 04 


9 

36602.20557 

0. 554 3 80621 £-07 

-0 .252246060E— 03 


10 

46692.62500 

0.12686081 0E-06 

-0 .9082 53 556 E- 01 


11 

51315. 99414 

0. 174881022E-06 

-0.101045303 


12 

56079. 86719 

0 • 2364 97 82 9E-06 

-0.10 2 4 74 4 7 3 


13 

60826. 1 8848 

0. 3117 392 1 0E-06 

-0 .9 9 7303 8 36E-0 1 


14 

63874.28467 

0. 368122151E-06 

-0. 105319239 


STATION NO 

Ml 

M2 

B1 82 


l 

500000000000 

500000000000 

-0 

0 

? 

0.855753 057 

191.1652775 

0.355753057 

191.1652775 

3 

0.8212 72 783 

185.7605033 _ 

3.378438413 

379.4372559 

4 

-208.6824379 

-47801 •14643 " 

-201.9848442 

-47416.86279 

5 

7,505106270 

1741. 419479 

-605.7670593 

-46264.91553 

6 

-42.4542 7990 

-9972. 145020 

-1041.049042 

- 56 796.81592 

7 

146. 4223499 

34758.61 182 

-1369.194733 

-22 706.71436 

8 

145. 896 0609 

34988.15723 

-1403.083771 

12021.40430 

9 

148. 4611 835 

35963.74561 

-1142.771103 

48119.230 47 

10 

37. 78408575 

9243.021851 

-698.6522064 

57890.10645 

11 

16. 5 6121421 

4096.619507 

-2 3 3 .9648 32 3 

62600.05176 

12 

17.54501176 

4385,690857 

320.6998253 

67632,77246 

13 

17. 50662494 

4421.687195 

376.0561447 

72736.79004 

14 

4. 635173559 

2350.302826 

1166 .672546 

76089.63965 
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The elastic-plastic program is quite short in running time. However, the creep pro- 
gram running time is dependent on the total number of time increments calculated. The 
example problem ran less than 1 minute for the elastic-plastic analysis and 2 minutes for 
the creep analysis. 
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APPENDIX B 


OPERATING INSTRUCTIONS 

The plastic flow and creep analyses were kept separate, for often only one analysis is 
used. Control cards were kept at a minimum. Each analysis will require new stress- 
strain and creep curves. The subroutine PLSTR, therefore, must be rewritten each time.. 
This subroutine occurs by the same name both in the elastic-plastic analysis and the creep 
analysis. The equations currently in PLSTR of the elastic -plastic solution are for the 
stress-strain curves of high-alloy steel at 1000°, 800°, and 600° F. The equations cur- 
rently in PLSTR of the creep analysis are the constant- temperature stress-strain - time 
relations of high- alloy steel at 1000°, 800°, and 600° F. 

Operating Instructions for "Stress" Code (Elastic- Plastic Analysis) 

The following is a summary of the input cards for the ’’Stress" code. They are 


numbered 1, 2a, 2b, etc. 

The letter designations are for continuation cards when the 

number of mesh points exceed a single-card limitation. 

Quantity 

Format 

Card columns 

Remarks 

Card 1 




NOFS 

15 

1-5 

Number of mesh points that the cylin- 
der is divided into 

NJ 

15 

6-10 

Control designation: 0 if cylinder con- 
tains a hole; 1 if cylinder is solid 

Card 2a, 

2b ... 



TEMP 

8F10.7 

1-10 

Temperature of each mesh point; addi- 



11-20 

tional cards where number of mesh 



21-30 

points exceeds 8, 16, etc. 



etc. 


Card 3a, 

3b ... . 



ALPHA 

8F10.7 

1-10 

Coefficient of thermal expansion cor- 



11-20 

responding to temperature at each 



21-30 

mesh point 



etc. 
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Quantity 


Format 


Card column 


Remarks 


Card 4a, 4b . . . 


E 

8F10. 7 

1-10 

11-20 

21-30 

etc. 

Modulus of elasticity corresponding 
to temperature of each mesh point 

Card 5 

XMU 

F10. 7 

1-10 

Poisson's ratio 

SIGRR 

F10. 7 

11-20 

Pressure load on outside radius of 
cylinder, psi 

SIGRI 

F10. 7 

21-30 

Pressure load on inside radius of cy- 
linder, psi; 0 if cylinder is solid 

Card 6 

NO FI 

15 

1-5 

Maximum number of iterations allowed 
(If NO FI is reached, program con- 
tinues to print output. If convergence 
is achieved prior to iteration limita- 
tion, program prints final solution. ) 

Card 7 


RAD F10. 7 1- 10 

11-20 

21-30 

etc. 


Radial location corresponding to each 
mesh point, in. 


Description of Output for "Stress" Code (Elastic- Plastic Analysis) 

The following is a brief summary of the output for the "Stress" code. Each Roman 
numeral refers to a new section of output. 
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Section 

Quantity 

Remarks 

I 

NO OF ITERATIONS 

Iteration limitation (NOFI) placed upon 
program (input) 


NO OF STATIONS 

Number of mesh points (NOFS) the cy- 
linder is divided into (input) 

II 

STATION NO 

Mesh point number 


CYL RADIUS 

Radial location (RAD) of mesh point 
(input) 


TEMPERATURE 

Temperature (TEMP) of cylinder at 
each mesh point (input) 


ALPHA 

Coefficient of linear expansion 
(ALPHA) of material at each mesh 
point (input) 


E 

Modulus of elasticity E of material 
at each mesh point (input) 

III 

POISSON’S RATIO 

Poisson’s ratio of material (input) 


RAD R STRESS 

Pressure at cylinder outer radius 
(input) (SIGRR) 


RAD I STRESS 

Pressure at cylinder inner radius 
(input) (SIGRI) 

IV 

STATION NO 

Mesh point number 


TOTAL STRAIN R 

Total radial strain at each mesh point 
(eq. (3» 


TOTAL STRAIN T 

Total tangential strain at each mesh 
point (eq. (4)) 


TOTAL STRAIN Z 

Total longitudinal strain at each mesh 
point (eq. (5)) 

V 

STATION NO 

Mesh point number 


PLASTIC STRAIN R 

Plastic radial strain at each mesh 
point (eq, (14)) 


PLASTIC STRAIN T 

Plastic tangential strain at each mesh 
point (eq. (15)) 


CYL RADIUS 

New location of mesh point due to 
radial strain 

VI 

STATION NO 

Mesh point number 


RADIAL STRESS 

Radial stress in cylinder at each 
mesh point (eq. (27)) 


TANGENTIAL STRESS 

Tangential stress on cylinder at each 
mesh point (eq. (27)) 
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Section 


Quantity 


Remarks 


VI 

LONG STRESS 

Longitudinal stress on cylinder at 
each mesh point 


EQUIV STRESS 

Effective stress at each mesh point 
(eq. (10)) 

VII 

(No label) 

Number of iterations at convergence; 
last calculated convergence tolerance, 
DELR, DELT, and constant C3 (See 
listing) 

vrn 

STATION NO 

Mesh point number 


ELASTIC STRAIN 

Equivalent elastic strain at each mesh 
point (eq. (12)) 


PLASTIC STRAIN 

Equivalent plastic strain at each mesh 
point (eq. (11)) 


P PRIME 

*Value of P^ of equation (22) at each 
mesh point 

IX 

STATION NO 

Mesh point number 


Ml 

* Value of M^ of equation (25) at 
each mesh point 


M2 

* Value of Mg n of equation (25) at 
each mesh point 


B1 

*Value of equation (30) at each mesh 
point for radial coefficients 


B2 

*Value of equation (30) at each mesh 
point for circumferential coefficients 

(*This information is printed out to aid in evaluating the validity of the solution. ) 


Operating Instructions for "Creep" Code (Creep Analysis) 

The following is a summary of the input cards for the "Creep” code. The cards are 
numbered similar to the stress code (i. e. , 1, 2a, 2b, etc. ). 
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Quantity 


Format 


Card columns 


Card 1 
NOFS 

15 

1-5 

NJ 

15 

6-10 

NI 

15 

11-1.5 

Card 2a, 2b . 

. . 


RADX 

8F10. 7 

1-10 



11-20 

21-30 

etc. 


Card 3a, 3b . . . 


TEMP 

8F10.7 

1-10 

11-20 

21-30 

etc. 

Card 4a, 4b . 

• * 


ALPHA 

8F10. 7 

1- 10 


11-20 

21-30 

etc. 


Card 5a, 5b . . . 

E 8F10.7 1-10 

11-20 

21-30 

etc. 


Remarks 


Number of mesh points that the cylin- 
der is divided into 

Control designation: 0 if cylinder con- 
tains a hole; 1 if cylinder is solid 
Control designation; printout of output 
will occur after NI time increments 


Radial location corresponding to each 
mesh point prior to any creep strain 


Temperature of each mesh point; ad- 
ditional cards where number of mesh 
points exceeds 8, 16, etc. 


Coefficient of thermal expansion cor- 
responding to temperature of each 
mesh point 


Modulus of elasticity corresponding 
to the temperature of each mesh point 
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Quantity 


Format 


Card columns 

Card 6a, 6b . . . 

EPR 8F10. 7 1-10 

11-20 

21-30 

etc. 

Card 7a, 7b . . . 

EPT 8F10.7 1-10 

11-20 

21-30 

etc. 

Card 8a, 8b . . . 

PP 8F10.7 1-10 

11-20 

21-30 

etc. 

Card 9 


XMU 

F10. 7 

1-10 

SIGRR 

F10. 7 

11-20 

SIGRI 

F10. 7 

21-30 

ACC 

F10. 7 

31-40 

Card 10 

NOFI 

15 

1-5 

NODT 

15 

6-10 


Remarks 


Plastic radial strain at each mesh 
point, obtained from the output of 
"Stress” code 


Plastic tangential strain at each mesh 
point, obtained from output of "Stress" 
code 


Value of equation (22) of "Stress" 
code at each mesh point, obtained 
from output of stress code (P PRIME) 


Poisson’s ratio 

Pressure on outside radius of cylin- 
der, psi 

Pressure on inside radius of cylinder, 
psi; 0 if cylinder is solid 
Acceleration factor for time incre- 
ment; initial time increment, 0. 01 hr 


Maximum number of iterations al- 
lowed (If NOFI is reached, program 
continues to print output. If converg- 
ence is achieved prior to iteration 
limitation, program prints final solu- 
tion. ) 

Total number of time increments for 
which creep calculations are performed 
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Description of Output for "Creep" Code (Creep Analysis) 


The following is a brief summary of the output for the ''Creep" code. Each Roman 
numeral refers to a new section of output. 

Section Quantity Remarks 


I NO OF ITERATIONS 
NO OF STATIONS 

II STATION NO 
CYL RADIUS 

TEMPERATURE 

ALPHA 

E 

III STATION NO 

PLASTIC STRAIN R 

PLASTIC STRAIN T 

P PRIME 

IV POISSON’S RATIO 

RAD R STRESS 

RAD I STRESS 

V STATION NO 

TOTAL STRAIN R 

TOTAL STRAIN T 


Iteration limitation placed upon pro- 
gram (NOFI in input) 

Number of mesh points cylinder is 
divided into (NOFS in input) 

Mesh point number 

Radial location of mesh points prior 

to strain (RADX in input) 

Temperature of cylinder at each mesh 
point (TEMP in input) 

Coefficient of linear expansion of ma- 
terial at each mesh point (ALPHA in 
input) 

Modulus of elasticity of material at 
each mesh point (E in input) 

Mesh point number 

Plastic radial strain at each mesh 

point (EPR of input) 

Plastic tangential strain at each mesh 
point (EPT of input) 

Value of equation (22) of "Stress” 
code at each mesh point (PP of input) 
Poisson's ratio of material (XMU of 
input) 

Pressure at cylinder outer radius 
(SIGRR of input) 

Pressure at cylinder inner radius 
(SIGRI of input), 0 if cylinder is solid 
Mesh point number 

Total radial strain at each mesh point 
(eq. (35)) 

Total tangential strain at each mesh 
point (eq. (36)) 
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Section Quantity 

V TOTAL STRAIN Z 

VI STATION NO 

INCH PLAS STR R 

INCR PLAS STR T 

CYL RADIUS 

STATION NO 
PLASTIC STRAIN R 
PLASTIC STRAIN T 
STATION NO 
RADIAL STRESS 

TANGENT STRESS 

LONG STRESS 

DC (No label) 

X STATION NO 

EQUIV STRESS 

PLASTIC STRAIN 

P PRIME 

XI STATION NO 
Ml 

M2 


Remarks 

Total longitudinal strain at each mesh 
point (eq. (37)) 

Mesh point number 

Plastic radial strain occurring during 
time increment AT at each mesh 
point (eq. (39)) 

Plastic tangential strain occurring 
during time increment AT, at each 
mesh point (eq. (40)) 

New radial location of mesh point due 

to radial strain 

Mesh point number 

Total plastic radial strain 

Total plastic tangential strain 

Mesh point number 

Radial stress on cylinder at each mesh 
point (eq. (27)) 

Tangential stress on cylinder at each 
mesh point (eq. (27)) 

Longitudinal stress on cylinder at each 
mesh point 

Number of iterations at convergence; 
last calculated convergence tolerance 
constant, C3, and total time T (See 
listing) 

Mesh point number 

Equivalent stress at each mesh point 

(eq. (10)) 

Equivalent strain due to creep at each 
mesh point (eq. (43)) 

*Value of equation (42) at each mesh 
point 

Mesh point number 
* Value of M 1 of equation (25) at 
each mesh point 

* Value of M 0 _ of equation (25) at 
n 

each mesh point 


VII 

vm 
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Section 

Quantity 

Remarks 


Bl 

* Value of equation (30) at each mesh 



point for radial coefficients 


B2 

*Value of equation (30) at each mesh 


point for circumferential coefficients 

(*This information is printed out to aid in evaluating the validity of the solution. ) 
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APPENDIX C 


PROGRAM LISTINGS 
Elastic- Plastic Program Listing 


$ I BJGB GO, MAP, SOURCE 

SI8FTC STRESS LIST, DEBUG 

DIMENSION C( 30) ,0(30} ,F i 30) ,G! 30 ) ,HT 130) , RAD (30 ) ,E ( 30 ) ,CP( 30), 
10PC30) ,FP ( 30 ) ,GP { 30) , HP (30) ,PP(30) ,RL1 1 (30 ) ,RL12 (30 I , RL21 < 30 ), 
2RL22I 30) ,RM1 ( 30) ,RM2(30) ,AR(30) ,AT(30) ,B1 ( 30 ) , B2 (30 ) , S IGR( 30), 

3SI GT( 30), S I G 2 ( 30) ,SISE (30) , TEMP (30) ,EPR(30) , EPT(30 ) ,ETR(30 ) , ETT( 30 
4) , EET ( 30) , EP ( 30) , ALPHA { 30) ,ETZ(30) 

COMMON TEMP, EET, EP f EPR ,EPT,RAD,I 

901 FORMAT ( 8F 10. 7) 

90 3 FORMAT (415) 

READ ( 5,903) NOFS,NJ 

READ (5,901) (TEMPI I ) ,i=l,NOFS) 

READ (5,901) (ALPHAU) ,1=1 ,NOFS) 

READ (5,901) < E( I 3,1=1 ,NUF S) 

READ (5,901) XMU , SI GRR , SI GR I 

902 READ ( 5,903) NUF I 

READ (5,901) (RADII ), 1=1, NQFS) 

904 FORMA T ( 20H1N 0 OF ITERATIONS = G13. 5 ,10 X , 1 7HN0 OF STATIONS = G13.5) 
WRITE (6,904) NOF I, NOF S 

906 FORMAT ( 1 1HK STATION NO , 1 5X , 1 1HC YL RAD I US, 9 X , 1 1 HTEMPERATURE, 9X , 
15HALPHA , 1 5X, 1HE ) 

WRITE (6,906) 

907 FORMA T ( 5G 20. 9) 

DO 908 I = 1 ,NGF S 

908 WRITE (6,907) I ,RAD( I ) , TEMP ( I) , ALPH All ) , E ( I ) 

910 FORMAT ( 18HJP 01 SSONS RATIO = G1 5 . 5 , 3 X ,1 5HRAD R STRESS = G15.5,5X, 
115HRAD I STRESS = G15.5) 

WRITE (6,910) XMU , SIGRR , SI GRI 

*************************************** 

BEGIN CALCULATIONS 

*********** *** *** ******************** ** 


SUMA =0.0 
DO 3 I =2, NOF S 
3 HT ( I)= RADII )-RAD ( 1—11 
HT ( 1) = H T ( 2 ) 

*************************************** 

CALCULATE THE GEOMETRICAL COEFFICIENTS 
*************************************** 


DO 1 J =1, NOF S 
IF ( J- 1 ) 2,99,2 
CALCULATE C 

2 C(J) = WHT4 J) + 1./(2.*RADU) ) 

CALCULATE D 

D( J ) = 1 • / ( 2 •♦RAD ( J) ) 

CALCULATE F 

FCJ) = i./HT(J)-l./(2.*RAD(J-l)) 
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on o o o o r> o no o o o o on on o o o o o o 


CALCULATE G 

GIJ) = 1 • / ( 2 .*R AD4 j-.l ) ) 

CALCULATE C PRIME 

CP I J) = - XMU/ IHTU) *E (J) j -XMU**2 / IH T I J I *EC J ) ) - ( 1 •+XMU ) / ( 2 • *EC J -V* 
IRAQI J ) ) 

CALCULATE Q PRIME 

DPI J) * l./IHTI J)^E(J)i-XMU«^2/CHT{J)«EU) ) + Cl .+XMU )/ l 2.*E(J ) *RAD 
1U)I 

CALCULATE F PRIME 

FP(J> * -XMU/IHT ( J ) *EU-1))-XMU**2/4HTU) *E I J~1 ) ) + C 1. +XMU) / 

1(2.*EU-1)*RADU-1)) 

CALCULATE G PRIME 

GPU) ~ 1./4HTIJ) *E U-l) )-XMU**2/ IHT4 J) *E U-i ) ) -( 1 . +XMU f / 

1( 2.*E4 J-l)*RAD(J-l)i 

CALCULATE H PRIME 

HPU) = -4 1. +XMU)/HTU ) #( ALPHA U)*T£MPU) -ALPHA U- 1 mEMPU-1) ) 

CALCULATE P PRIME 
PP(J) = 0.0 

DEBUG C(J )iO{J)i Ft JltGfJI 

DEBUG CPU) ,DPU) ,GPI J) ,FPI J) r HPU) 

CALCULATE MATRIX COEFFICIENTS LI1,L12 
DEN = (CIJ1*0P{J)+CPIJ)^{J)) 

R L 1 1(J ) = 1 DP U ) *F (J ) + D ( J J *F P U ) ) /D E N 
RL12IJ) ^4QPU)*GU)+DU>*GPU) ) /DEN 
RL 2 IU I =4 Cl Jj*FP(J I-CP U) *FU) I /DEN 
RL 221 J ) *4 C U ) *GP ( J )-CP I J) *G (J)) /DEN 
RM IU ) =D(Ji *4 HP U )+PP U ) ) /DEN 
RM2U) <U)*lHPUHPP(J))/OEN 
DEBUG RL1 14 J ) f RLL2 ( J) ,RL21U) ,RL22U> 

CALCULATE COEFFICIENTS AR , AT 

AR4J ) =RL 114 J ) *AR l J-l J *RL12 ( J) *A TOl) 

ATI J) =RL21I J)*ARU-1 ) +RL22 I J) *ATU-1 ) 

DEBUG ARC if, AT ( J ) 

GO TO 1 

9 9 IF4NJ) 6,6,7 

SOL ID CYL INOBR 
1 ARC J ) = 1.0 
ATI J ) - 1.0 
Bl(J) = 0.0 
B24 J ) = 0.0 
GO TU 1 

CYLINDER WITH HOLE 
6 ARU) ^ 0.0 
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ATCJI = 1.0 
BHJ) = -SIGRI 
82(J> = 0.0 
C 

1 CONTINUE 
C 

61 DO 63 IP =1, NOFS 
EPUP) = 0.0 
EPR ( IP ) = 0.0 
EPT(IP) = 0.0 
ETR(IP) = 0.0 
ETT(IP) = 0.0 
63 CONTINUE 
C 

C***************** ********************** 

C ITERATION LOOP 

£*****************#*******************& 

DO 501 NIT= 1 , NOF I 
SUMA = SUMA + 1.0 
NK « 0 
XC4 = 0.0 
XC3 = 0.0 
504 DO 9 K =2 » NOF S 
C 

DEN= ( Cl K I *D P ( K J *C P ( K ) *D ( K ) ) 

RM 1 1 K ) = OIK) * (HP ( K ) +P P ( K > I /DEN 
RM 2 ( K I = Cl K) * (HP {.KJ+P.P (K) ) /DEN 
C 

c 

10 B11K|= RL I K K)*B1( K-l ) +RL12 ( K) *B2(K-1)+RM1 (K ) 
82 ( K 1= RL2K KI*B1 (K-1 »+RL 22 (K) *B2< K-1 ) *RM2 CKI 
15 CONTINUE 
9 CONTINUE 
C 

DO 39 1= 1 « N OF S 
C 

C CALCULATE RADIAL AND TANGENTIAL STRESSES 
C 

KPT= NOFS 
IF( 1-1) 20,25,20 
C 

25 IF (NJi 24,24,23 
C 

C CENTER STRESSES FOR A SOLID CYL 
C 

23 S I GT( II = -B 1 (KPT) /AR( KPT) -SIGRR/ARI KPT) 
SIGR(l) = SIGT(l) 

GO TO 29 
C 

C STRESSES FOR HOLED CYL 

C 

24 SIGRC 1 )=- SIGRI 

SIGH!)* -B1 (KPT) /AR( KPT I - SIGRR/ ARIKPT I 

GO TO 29 
C 

20 IFCI-KPTI 21,22,21 
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22 SIGRIKPT)=-SIGRR 
GO TO 28 
C 

21 SIGRC I )= ARC I)*$IGTC 1)+B1 II) 

28 SIGT ( I ) = AT{ I)*$IGT<1)+82(II 
C 

29 CONTINUE 
C 

S I GZ ( I ) * XMU^CSIGTCn + SIGRCn I - EUr*ALPHAUI*TEMPm+E{I J*EPRU 
1)*£4 f )*EPTI I I 
39 CONTINUE 

DO 35 I - 2, NOF S 

XC2 = (EC I)+E( I-1>)/4.*CRADCI)**2-RADSI-1>**2> 

XC4 * XC4+XC2 

XC1 * (SIGZt I J + S IGZ ( I - 1 1) /4. *C RAO m**2- RAO (I-i)**2l 
35 XC3 * XC3 + XC1 

C3 = CXC3 - ( SIGRI * RADII) **2 ) /2 • 0 ) /XC4 
DO 37 1-1 ,NGFS 

37 SIGZin * SIGZU) - C3 *E(l) 

DO 30 1=1 » NOF S 

C CALCULATE THE TOTAL STRAINS 
C 

ETRX= ETR I £ ) 

ETTX= ETT ( I ) 

C 

ETRI I ) = C 1. /EC U 1*1 SIGRC I )-XMU* (SIGT ( I) + S IGZ ( I ) ) ) + ALPHA! I ) *TEMP( I 
I ) +EPR ( I) 

ETT(I) = Cl. /EUn*tSlGT<l)~XMU* (SIGRC I) + S IGZ (I) > > +ALPHAI I)*TEMP( I 
1I+EPT1 II 

£ T Z ( I ) = X I. /EC II J*C SIGZM1-XMU* 4 SlGRll I +SIGT! I H)+ALPHA{ U*TEMP 
II I l-EPRC I J-EPT4 1 ) 

C 

C CALCULATE THE EQUIVALENT STRESS 

PRA- C (SIGRC I )-SIGTCI I 1**2+ {SIGRC It-SIGZ1I )> *«2+(S IGT C U-S IGZC I ) ) 
1 ** 2 ) 

S IGEC I )= ( 1 • / 1 • 41 42 1 ) * SQR T C P RA I 
C 

C* ********** **************************** 

C TEST FOR CONVERGENCE 

C * ************************************** 

DELR- CETRX-ETRM n/ETRCU 
DELT= CETTX-ETTCin/ETTCn 
C 

IF (A8SC0ELR )~.001) 31*31,32 

32 NK= 1 

31 IF (ABSCDELT)-.OOl) 33,33,34 
34 NK - 1 

33 CONTINUE 
C 

C CALCULATE THE EQUIVALENT TOTAL STRAINS 

STA = CCETRC I ) -ETT Cl ) t **2* CETT C I T-ETZC 1 1)**2* < ETZC n-ETRC I) 1**2) 
■EE-Tli )» C 1. 41421/3. G)*SQRT(STA) 

C 

C 

C 

CALL PLSTR 
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CALCULATE PLASTIC STRAINS 

66 EPR{ I T= EP(I )/(3.*EET( I) )*(2.*ETRI I)-ETT( I )-ETZ{ I) ) 

EPT ( I ) = EP ( I ) /( 3.*EET( I) )*(2.*ETT(I) -ETR( I l-ETZC I ) ) 

500 CONTINUE 
IF(I-l) 71,3 0,71 

CALCULATE NEW VALUE FOR P PRIME 

71 PPIII ■« l-l./HT( I)-1./(2.*RAD( I) ) + XMU/HT(I))*EPT(IH-(l./(2.*RAD( I) 
1) + 

1XMU/HTCI> l*EPR(l)+(l. /HT 1 1 ) — 1. / ( 2. *RAD (1-1 ) T-XMU/HT 1 1) )*EPT( I-l) + 
21 l./(2.*RAD< I-1J I-XMU/HT ( I ) )*E PRII-1I 

30 CONTINUE 

IF(NK) 50 1 , 5 C3 , 5 01 

501 CONTINUE 
503 CONTINUE 

00 75 I = 1, NOF S 

75 RADI 1 1 = RAO ( I ) *11. + E TT (1)1 
C 

Q*4*************** ********************** 

C BEGIN PRINT OUT OF OUTPUT 

C* ** ******** ****** ********************** 

912 FORMAT (11HLSTATI0N NO , 1GX , 14HT0TAL STRAIN R, 6X , 14HT OTAL STRAIN T,6 
IX, 14HT0TAL STRAIN Z) 

WRITE 16,912) 

914 FORMAT (4G2Q.9) 

DO 913 I=1,N0FS 

913 WRITE (6,914) I ,ETR( I ) ,ETT( I ) ,ETZ( I } 

915 FORMAT ( 11HK STATION NO ,10X,16HPLASTIC STRAIN R, 4X, 16HPLAST I C STRAI 
IN T , 10X, 1 OHC YL RADIUS) 

WRITE (6,915) 

919 FORMAT (5G20.9) 

DO 918 I=1,N0FS 

918 WRITE (6,919) I ,EPR( I ) ,EPT ( I ) ,RAD( I ) 

920 FORMAT ( 11HKSTA TION NO , 10X .13HRADIAL STRESS ,7X , 14HT ANGENT STRESS, 6X 
1, 11HL0NG STRESS, 9X.12HEQUIV STRESS) 

WRITE (6,920) 

00 925 1=1, NOF S 

925 WRITE (6,919) I , SI GR( I ) , SIGT ( I ) , SIGZ ( I ) ,S IGE ( I ) 

WRITE (6,919) SUMA,DELR,DELT,C3 

929 FORMAT (11HKSTATIQN NO ,10X ,14HELASTI C STRAIN, 6X, 14HPLAST IC STRAIN, 
16X, 7HP PRIME ) 

WRITE (6,929) 

931 FORMAT (5G20.9) 

DO 930 1=1, NOF S 

930 WRITE (6,931) I ,EET( I ) ,EP( I ) ,PP( I) 

928 FORMAT ( 1 1HKSTATI ON NO , 10X ,2HM1 , 18X,2HM2 , 18X, 2HB1 , 18X , 2HB2 ) 

WRITE (6,928) 

DO 936 1 = 1, NOF S 

936 WRITE (6,931) I ,RM1( I ) ,RM2 (I) ,B1(I),B2(I) 

GO TO 902 
END 
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SIBFTC PLSTR LIST, DEBUG 

SUBROUTINE PLSTR 

DIMENSION Cl 30) ,0130) »FI 30) ,G( 30 1 ,HT <301 ,RAD (3G ) ,E 1 39 > ,CP( 30), 

IDP (301»FP(30) ,GP ( 30) ,HP( 30) ,PP(30) ,RL11(30 ) , RL12 (30 ) , RL2K 30), 
2RL22( 30) ,RMi (30) ,RM2(3Q) ,AR(30) »AT(30I ,81 (30 ) , B2 ( 30 ) ,S IGRI30), 

3SI GT ( 30), SIGZI 30) ,SIGE (30) , TEMP (30) , EPR(30 ) , EPT ( 30 ) , ET R( 30 ) , ET T ( 30 
4),£ET( 30) ,EP (30) , ALPHA (30) ,ETZ(3GI 
COMMON TEMP, EET,EP,EPR ,EPT , RAD-,1 
C 

C******** ****** ************************* 

C MATERIAL IDENTIFICATION 

£*** *********** *** ********************** 

C 

IF (RADUJ - RAD ( 14) } 10,10,11 
11 GO TO 500 
C 

£* ** ************** ********************** 

C TEMP IDENTIFICATION MAT 1 

£*** *********** ************************* 

C 

10 IF (TEMPI II - 1000.) 25,20,20 

25 IF (TEMP(I) - 800. I 35,30,30 
35 IF ( TEMPI I 1 - 600.) 45,40,40 
45 IF I TEMPI 1 1 - 400.) 55,50,50 
55 GO TO 50 

C 

C STRAIN CALC TEMP 1000 

C 

20 IF (EET(I) - .0045) 26,26,27 

26 EP(l) = • 6*E£T( I 1 - .0012 
GO TO 400 

C 

27 EP ( I ) = .98*EET(1) - .0027 
GO TO 400 

C 

C STRAIN CALC TEMP 800 

C 

30 EPU) = 0 . 922 * EET(I) - .00207 
GO TO 400 
C 

C STRAIN CALC TEMP 600 

C 

40 EPU) * .898 * EETII) - .00184 
GO TO 400 
C 

C STRAIN CALC TEMP 400 

C 

50 IF (EETII) - .0017) 51,51,52 

51 EPU) * .81 * EETII ) - .00154 
GO TO 400 

C 

52 EPII) = .931 * EETII) - .00202 
400 IF (EPII ) .LT. 0.0) GO TO 80 

GO TO 500 
C 

c 

80 EPU) * 0.0 
EPRI I ) * 0.0 
EPT (I ) = 0.0 
500 RETURN 
END 
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Creep Program Listing 


*IBJOB GO, MAP, SOURCE 

SlBFIC CREEP LIST, DEBUG 

DIMENS ION C( 30) ,0(30) ,F(30) ,G(30),HT(30) ,RAD ( 30 ) ,E < 30 ) ,CP< 30), 

1DPI 30) ,FP ( 30) ,GP ( 30) ,HP(30) ,PP(30) ,RL1 1 (30 ) , RL12 ( 30 ) , RL21( 30 ), 
2RL22( 30 ,KM1 ( 30 ) ,RM2(30) ,AR(30) ,AT (30) ,BU3U ) , B2 ( 30 ) ,S IGR( 30), 

3S I GT( 30),SIGZ(30) , SI GE (30) , TEMP (30) ,EPR( 30 ) , EPT (30 ) , ET R( 30 ) , ET T ( 30 
4),ALPHA(30),ET2(30) ,0EPR(30) ,OEPT(30) ,OEP (30 ) , RADX { 30 ) 

COMMON TEMP, DEP, SIDE, I ,DT,kAO 

901 FORMAT (8f 10. 7) 

90 3 FORMAT (415) 

READ 15,903) NOF$,NJ,NI 

READ (5,901) (RADX(I) ,I=i,NQFS) 

READ (5,901) (TEMPd ) ,1=1 ,NOFS) 

REA u (5,901) (ALPHAU ) ,1=1 ,NOFS) 

READ ( 5, 90 i. ) (E(l) ,I = 1,N0FS) 

READ ( 5,901) (EPR(I), 1=1 ,NOFS) 

READ (5,901) (EPT(I), I=1,NUFS) 

READ (5,901) (PP(I), 1=1, NUF S) 

READ (5,901) XMU, S1GRR , SI GRI , ACC 

902 READ (5,903) NUFI.NUDT 

904 FORMAT (20H1NU UF iTEKATIUNS = G13.5 ,10X,17HNU OF STATIONS = G13..5) 
WRITE (0,904) NUFI.NOFS 

906 FURMAT (11HRSTATIUN NO ,1 5X ,11HC YL RAD I US , 9X , 1 1HTEMPER ATUR E, 9X , 
15HALPHA, 15X, 1HE) 

WRITE (0,906) 

90 7 FORMA T ( 5G20. 9) 

DU 908 1 = 1, NUF S 

908 WRITE (6,907) I ,RADX( I ) , TEMP ( I ) , ALPHA ( I ) , E ( I ) 

980 FORMAT { 11HK STATION NO ,10X ,16H PLASH C STRAIN R ,4X, 16HPLAST IC STRAI 
IN T , 4X , 7HP PRIME) 

WRITE (6, 980) 

981 FORMAT (5G2U.9) 

DO 982 I = l.NUFS 

982 WRITE (6,981) I.EPR(I), EPT(I), PP ( I ) 

910 FORMAT U8HJP 01 SSONS RATIO = Gi 5 . 5 ,3X ,1 5HRA0 R STRESS = G15.5.5X, 
115HRA0 I SIRESS = G15.5) 

WRITE (6,910) XMU,SIGRR,SI GRI 

*#*** 4 =**#*#***#***** ****** 

BEGIN CALCULATIONS 

DT = .01 
T= 0.0 
NH = NI-1 

61 DO 63 IP =1, NOFS 
ETR( IP ) = 0.0 
£TT( IP ) = 0. 0 
DEP (IP) = 0.0 
DLPR( IP)= 0.0 
DEPTI IP)= 0.0 
63 CONTINUE 
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C iIMfc INCREMENTS LOUP 

J, *** *** ****** ******************** 

C 

DO 505 J J =i» NUDT 
SUMA =0.0 
DT= DT*ACC 
WH * NH+1 

6 7 DO 68 IJ = 1 fNOFS 

£PRU0) = DEPRIIJ)^ EPRil J J 
tpruj)= OEP T ( UJ+EPTUJ) 

c 

68 CONTINLE 

DU 75 I=i,N0F5 

75 RADII ) = RAD X( I I* 4 1 • + ETTUii 
DU 3 l *4 9 NOFS 
3 HT4 i)« RAD 4 I 1-RAD4 1-1 \ 

Hlii> = HU2) 

****** *** *********************** 

CALCULATE THE GEUMETRI CAL COEFFICIENT'S 

Xf^i:****#*****#*****^************ 

DU I J=1,NUFS 
IF(J-I) 2,99,2 
CALCULATE C 

2 CU) = i./HTUVU./I2.*RA0iJJ ) 

CALCULATE D 
DU) = I. / 1 2 »*RAD (J) ) 

CALCULATE F 

H J I = I./HT < J)-1./I2.*RAD( J-l) ) 

CALCULATE G 

GU > * l./(2.*RA0U-i) ) 

CALCULATE C PRIME 

CPU) * -XMU/(HTU )*EU) )-XMU**2/<rtTU)*EU) > -I1« + XMU)/{2.*EU )* 
1RADU ) ) 

CALCULATE D PRIME 

DP ( J ) = l./(HTU)*£U))-XMU**2/ (HT IJ)*EU) l+tl.+XMU )/( 2.*E(«J )*RAD 
11JI) 

CALCULATE F PRIME 

FPU) * -XMU/1HTU) *EU-l) )-XMU**2/4HTU) *E U-l ) ) + ( l . + XMU) / 

It 2.*E( J-l )*RADU-l) ) 

CALCULATE G PRIME 

GPU) = 1./IHTU) *EU-1) )-XMU**2/(HT (J) *E{ J-l) ) -( 1 ,+XMU )/ 

l(2.*Et J-i)*RADU-l)) 

C 

; CALCULATE H PRIME 

HP (J ) * - I l»+XMU)/HTU)*t ALPHA t J ) *T£MP IJ) - ALPHA U-l ) #TEMP( J-i) ) 

C 

DEBUG CU ),0U) ,P(J) ,GU) 
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DEBUG CP ( J ) , DP 4 J ) , GP U J * F P 4 J ) ,HPU) 

C CALCUL ATE MATRIX COEFFICIENTS Lil,L12 

DEN * 4C4 JJ*DP4 Jl+CPC JJ*D( J)I 
Rtili J i * I DP < J J *F < J ) +D <J ) *FP 141 ) /DEN 
R L 1 2 4 J ) = 4DP U)*G4 JL+04u)*GP4 JI.I/DEN 
RL214 J ) ~4C( J)*FP4 JJ-CP4 J)*F4.JJ ) /DEN 
Rc 224 J i = 4C( J)*GP4JHCP( JJ*G4J) ) /DEN 
RM 14 J > =04 J)*(HP-(-J)+PP4 J) ) /DEN 
RM24 J } =C(J)«( HP 4 J ) +P P 4 J 1 J /DEN 
DEBUG RCii( J J »RL12 4 J ) ,RL2i ( J) ,RL22 4 J) 

C 

C CALCULATE COEFFICIENTS AR » AT 

AK4 J ) =RLlLi JJ*AR4 j-l)+RL124 J)*AT4 J-l ) 

A T 4 J ) -RL21C Ji *AR4 U-l ) +RL2 2 ( J > *A T 4 J~ 1 ) 

DEBUG AR4 Jit AT4J) 

C 

GO TO 1 
C 

95 1F4 NJ ) 6» 6 * 7 
C 

C SOLID CYL 

C 

7 AR 4 J i ^ 1.0 
AT 4 U I * i.O 
B14 J ) = 0.0 
B 2 4 J ) = 0.0 
GO TO 1 
C 

C HOLED CYL 

C 

6 AR4J) = 0.0 
AT 4 J 1 - 1.0 

= -SIGRI 
B24 J) = 0.0 
C 

1 CONTINUE 
C 

C **##$:**** 

C ITERATION LOOP 

C ******************************** 

c 

66 DO 501 NIT= 1,NUFI 
SUM A = SUMA + 1.0 
NK = 0 
XC4 '0.0 
XC 3 = C.O 
504 Ou 9 K=2» NUF S 
C 

OEN= ( C(K J*DPIK)+CP(K>*0(K) ) 

KHKK i = 0 ( K J #( HP (Kl +PP tK) ) /OEN 
Rrt2(K ) = C(M*(HP(K)+PP(K))/OEN 

C 

IQ B 1 ( K } = RL 1 1 ( K ) *B 1 { K- 1 J+RLI ? ( Kl *.8.2 ( K -1 ) +RMI ( K ) 
8 2 ( K ) = RL 21 ( K ) *B 1 < K- 1 ) +RL2 2 1 K ) *B 2 1 K- 1) +R M2 ( K > 
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15 CONTINUE 
9 CONTINUE 

00 09 1= l.NUFS 

CALCULATE RADIAL AND TANGENTIAL STRESSES 


KP T= NOES 

IF I 1-1 » 20,25,20 

25 IF CNJ) 24,24,23 


CENTER STRESSES FUR A SOU 0 CYL 

23 SIGTIi) * -GH KPT ) /ARIKPT) -SIGRR/ARIKPTJ 
S luKI 1 I = Si GT 1 1 ) 

GO TO 29 


STRESSES FUR HOLED CYL 


24 SIGRl 1 )=- SIGR1 

Si GTI 11= — B1 (KPT) /AR I KPT) - S IGRR/ AR( KPT ) 
GO TO , 29 

20 IF I I-KPT) 21 ,22 ,21 

22 SIGRIKPT) =-SlGRR 
GO TO 23 

21 SIGRl I I* ARUI*SIGTI U+aili) 

23 SIGTi I )« ATI I ) * S I G T 1 1 ) +8 2 1 1 ) 


2 9 CONTINUE 

SIG2II) = XM U*( S I G T 1 1 ) + S I G R 1 1 ) ) - E 1 1 ) #ALPHA( I ) *TEMP( I ) *fc.{ I ) *EPR( I 
1J+EII )*EPTII I 
39 CONTINUE 

DU 35 I * 2, NUF S 

XC2 = IE I !)•»-£( i-l) ) /4. * IRAQ 1 1 ) **2-RA0( I-l ) *#2 ) 

XC4 = XC4+-XC 2 

XC1 = ISIGZII) + S IGZ l I-l) ) /4. * I RAD 1 1 ) 4*2— RAD I I -I ) **2 ; > 

3 5 XC 3 = XC3 «• XC1 

C3 = I XC 3 - ISIGRI * RAD III **2) /2 . 0 ) / XC4 
DO 37 1=1 ,NUF S 

37 S IGZ I I I * SIGZII) - C3 *EU> 

DU 30 1= 1 ,NOFS 

CALCUCATt THE TOTAL STRAINS 

ETRII ) = 1 1. /EIU J* (Si GR( I )-XMU* (SIGT ( I ) +S1GZ l ID) +ALPHAII) *TEMP( I 
ll+EPRI I) + DEPRin 

EITII ) = Il./EII) )* (SIGTII i-XMU* ISIGRI II +S IGZ I I))) + ALPHA! I )*TEMP( I 
1H-EPTI II *-DEPT 1 1 ) 

£12(11* I I ./El m*I SIGZI I )-XMU* I SIGR(I H-SIGTIH I I + ALPHA! 1 1 #T EMP 
1U I-EPRI I I-EPTI I J -DEPTH ) - DEPKII) 

LALCUL ATE THE EQUIVALENT STRESS 

PRA* i (SlGRI 1 1 -SIGT 1 1 ) )**2+ (SI GR 1 1) -SI GZ 1 1 ) i **2 + IS IGT I U-S IGZ 1 1 I I 
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i** 2 > 

S l GEL I )= I 1* /l *41421) #SQRT C PRA) 
C 

CALL PLSTft 


CALCULATE PLASTIC STRAIN INCREMENTS 
EPKX = DfcPRI I) 

EPTX = UfcPTUJ 

OEPR< I )= OEP ( i ) ZI2.fr SIGE (I) ) *1 2. *S ICR { I) -S IGT i 11 -S IGZC I) ) 
DLP T i i )= UEP m/(2.*SIG6( I ) ) * i 2. *S IGT ( U -S IGR< Il-S IGZl I I ) 

####### 

TEST FCR CUN VERGE NCE 
***#***##£#£ 

□ ELK = IEPRX-DEPRU) )/DEPRCI) 

OELT = 1EPTX-DEPT1 I ) )/DEPT I I ) 

5U0 CONTINUE 

IF tAiiStUELR J-.UOl) 31,31,32 
32 NK = 1 

31 IF IABSID£LT)-.001) 33,33,34 
34 NK = 1 
3 3 CUN TIN UE 


IF! 1- 1 ) 71,30, 71 

CALCULATE new value for p prime 

71 PPCI)= 1 - 1 ./BT 1 1 )- 1 • / 1 2, * RADI I ) ) +XMU/HT II) ) *( EPT (I ) +OEPT ( 1)1 + 

11 i./(2.*RA01 I) ) +XMU/H Til) ) * 1 EPR II) +DEPR1 I ) )+ l l./HT 1 I)-1./12.*RAD 
21 1-1) )-XMU/H T1 I) )# 1EP T (I -1 ) +DEPT II -1 ) ) + 11. /1 2 . #RAD1 I-lU-XMU/ 
3HT1 i) ) * 1 E P R 1 I— 1 J +DEPK1 1-1) ) 


30 CONTINUE 

IF INK) 501,5 03,501 
501 CONTINUE 
5u 3 Cun t i n uE 
T = T+OT 

IF1NI.NE.NH) 00 TO 505 
NH = 0 

****** ****** ********* *********** 

BEGIN PRINT OUT OF OUTPUT 

****** ****** ******************** 

912 FORMAT (ilHLSTAIIUN NO , 10X , 14HT OT AL STRAIN R,6X, 14HT0TAL STRAIN T,6 
IX, 14 HTUTAl STRAIN 2) 

WRITE 16,912) 

914 FORMAT 1*020.9) 

DU 913 I = 1 , N OF S 

913 WRITE 16,914) I ,E TR 1 1 ) ,E TT 1 1 ) ,E TZt I ) 

915 FORMAT 1 ilHK STATION NO ,10X , 16HI NCR PLAS STR R,4X,16HINCR PLAS STR 
1 T , 4X , luHCYL RADIUS) 
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WRITE 16,915) 

919 FORMAT (5G2Q.9) 

00 913 i = l,NGFS 

918 WRITE (6,919) I ,OEPRU ), DEPT (I) , RADII) 

960 FORMAT 111 HR STATION NO ,10X,16H PLASTIC STRAIN R,4X, .16 HP LAST (C STRAI 
IN T) 

WRITE (6,960) 

961 FORMAT (5620.9) 

00 962 1= 1, HUES 

962 WRITE (6,9oi) I ,EPR ( 1 ) ,EPT ( I ) 

92 C FORMAT (HHRSTATiON NO , 10X » 13HRA0IAL STRESS ,7X , 14-HT ANGENT STRESS, 6X 
1, 11HL0NG STRESS) 

WRITE (6, 920) 

00 925 1= 1 ,NUF S 

925 WRITE (6,919) I , S TOR < 1 ) , SI GT II) , SI GZ 1 1) 

WRITE 16,919) SOMA ,DELR ,OELT ,C3 , T 

929 FORMAT ( 11HR STATION NO ,10X , 14HEUUIV. STRESS , 6 X , 14HPLAST IC STRAIN, 
16X , 7HP PR I ME ) 

WRITE 16,929) 

931 FORMAT (5G20.9) 

DU 930 1=1, NOP S 

936 WRITE 16,931) I , SIGEII I ,0E P( I ) ,PP1 I) 

928 FORMAT ( 1 1HKSTA TION NO ,10X ,2HM1 , 1.8 X ,2HM2 , 18X , 2H81 , 18X , 2HB2 ) 

WRITE (6,928) 

OU 936 I = 1 ,NUF S 

936 WRITE (6,931) I ,RMi(I) ,RM2 M ) ,61 (I ) ,82 (I) 

50 5 CONTINUE 
GO TO 902 
END 
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$ I BF TC PLSTR LIST, DEBUG 


SUBROUTINE PLSTR 

DIMENSION Cl 30) ,01 30) ,F(30) ,G ( 30 ) , HT (30 ) , RADI 30) ,E (30 ) ,CP( 30), 

10P ( 30 ) ,FP ( 30 ) , G P ( 3 0 1 , H P < 3 0 ) , PP ( 3 0 ) , R LI 1 (30 ) , RL12 ( 30 ) , RL21 ( 30 ) , 
2RL22I 3 0 , RM1 (30) ,RM2( 30) ,AR(30) ,AT (30) ,B1 (30 ) , B2 ( 30 ) ,S IGR( 30 ) , 
3SIOf( 30) , SIGZ( 30) , SIDE (30) , TEMP (30) , EPR( 30 ) , EPT ( 30 ) , ETRI 30 ) , EfT ( 30 
4 1 , AL P H A I 3 0 ) , E T l ( 3 0 ) , 0 £ PR (30) ,DEPT(30) ,DEP( 30 ) , R AUX ( 30 ) 

COMMON TEMP, OEP, SIDE, I ,0T , RAD 

****** *** *********************** 

MATERIAL IDENTIFICATION 
****** ************************** 

IF (RAD(I) - RAO (14) I 10,10,11 
11 GU TO 500 

****** *** *********************** 

TfcMP i UEN J IF I.C ATI ON MAT I 

10 IF (TEMPI I ) - 1000.) 25,20,20 
25 IF (TEMPI I ) - 800.) 35,30,30 
35 IF ( TEMPI I ) - 600.) 45 ,40,40 
45 GO TO 40 

STRAIN CALC TEMP 1000 

20 OEP ( I ) = 4.4QE-17*SIGEII )**2.57 *DT 
GO TO 400 

STRAIN CALC TEMP 800 
30 OEP ( 1 ) = 8.0E-2i*SIG£(l)**3.Gl *DT 
GO TO 400 

STRAIN CALC TEMP 600 
40 DEP(I) = 1.0E-24*SIGE(I>**3.4 *DT 

400 CONTINUE 
50 C RETURN 
END 


46 



REFERENCES 


1. Millenson, M. B. ; and Manson, S. S. : Determination of Stresses in Gas-Turbine 

Disks Subjected to Plastic Flow and Creep. NACA TN 1636, 1948. 

2. Johnson, Donald F. : Analysis of Elastic- Plastic Stress Distribution in Thin-Wall 

Cylinders and Spheres Subjected to Internal Pressure and Nuclear Radiation Heating. 
NASA TND- 271, 1960. 

3. Mendelson, A. ; Hirschberg, M. H. ; and Manson, S. S. : A General Approach to the 

Practical Solution of Creep Problems. J. Basic Eng. , vol. 81, no. 4, Dec. 1959, 
pp. 585-598. 

4. Mendelson, Alexander: Plasticity; Theory and Application. Macmillan Co. , 1968, 

pp. 119-121. 

5. Conway, J. B. : Properties of Some Refractory Metals. V: Numerical Methods for 

Creep and Rupture Analyses, Rep. GEMP-397, General Electric Co. , Jan. 31, 1966. 

6. Dorn, J. E.; and Shepard, L. A. : What We Need to Know about Creep. Symposium 

on Effect of Cyclic Heating and Stressing on Metals at Elevated Temperatures. Spec. 
Tech. Publ. No. 165, ASTM, 1954, pp. 3-30, 

7. Yalch, J. P.; and McConnelee, J. E. : Plane Strain Creep and Plastic Deformation 

Analysis of a Composite Tube. Nucl. Eng. Des., vol. 5, no. 1, Jan. /Feb. 1967, 
pp. 52-62. 


E-4488 


NASA -Langley, 1969 32 


47 





